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We present a new measurement of the e + e~ — > J/ip cc cross section where the cc pair can fragment 
either into charmed hadrons or a charmonium state. In the former case the J/ip and a charmed 
hadron are reconstructed, while the latter process is measured using the recoil mass technique, which 
allows the identification of two-body final states without reconstruction of one of the charmonia. 
The measured e + e~ — > J/ipcc cross section is (0.74 ± 0.08 Ip'og) pb, and the e + e~ — > J/ip X non . C c 
cross section is (0.43 ± 0.09 ± 0.09) pb. We note that the measured cross sections are obtained from 
a data sample with the multiplicity of charged tracks in the event larger than four; corrections for 
the effect of this requirement are not performed as this cannot be done in a model-independent way. 
The analysis is based on a data sample with an integrated luminosity of 673 fb -1 recorded near the 
T(4S) resonance with the Belle detector at the KEKB e + e~ asymmetric-energy collider. 

PACS numbers: 13.66.Bc,12.38.Bx,14.40.Gx 



Prompt charmonium production in e + e~ annihilation is important for studying the interplay between perturba- 
tive QCD and non-perturbative effects. The production rate and kinematic characteristics of J/ip mesons in e + e~ 
annihilation are poorly described by theory, and even the production mechanisms are not understood. An effective 
field theory, non-relativistic QCD (NRQCD), predicts that prompt J/ip production at ^/s~10.6 GeV is dominated 
by e + e~ — * J/ip gg with a lpb cross section [lj; the e + e~ — * J/ip g contribution, which may be of the same order, 
is uncertain due to poorly-constrained color-octet matrix elements |2|. The e + e~ — > J/ipcc cross section is predicted 
to be ~ 0.05 — 0.1 pb 3], only ~ 10% of that for J/ipgg @]. (The estimate of the ratio is more precise, as QCD 
uncertainties partially cancel.) By contrast, Belle observed the ratio of the J/ip cc and inclusive J/ip production cross 
sections to be 0.59^^3 ± 0.12 [5J, and thus found <r(e + e~ — * J/tpcc)/a(e + e~ — > J/ipgg) > 1. Such a large value 
cannot be explained within the NRQCD framework, however some alternative approaches (see e.g. Ref. |6|) can 
accommodate it. 

In this report we present a new measurement of the e + e~ — > J/ip cc cross section. This process can be experimentally 
tagged by the presence of another charmed particle (either charmonium or charmed hadrons) in the event in addition 
to the reconstructed J/ip. The technique used in this analysis allows the model dependence of the result to be 
removed, reducing the systematic uncertainties. Production of the J/ip via mechanisms other than e + e~ — ► J/ipcc is 
also studied. The J/ip momentum spectrum, and helicity and production angle distributions, are measured for both 
e + e~ — * J/ipcc and J/ip X non _ cg processes. The analysis is performed using data recorded at the T(45) and in the 
continuum 60 MeV below the resonance, corresponding to integrated luminosities of 605 fb -1 and 68 fb^ 1 , respectively. 
The data are collected with the Belle detector 0] at the KEKB asymmetric-energy e + e~ collider 0]. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector (SVD), a 50- 
layer central drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), a barrel- like arrangement 
of time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located outside 
the coil is instrumented to detect mesons and to identify muons (KLM). Two inner detector configurations were 
used. A 2.0 cm beampipe and a 3-layer silicon vertex detector were used for the first sample of ~ 156 fb -1 , while a 
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1.5 cm beampipe, a 4- layer silicon detector and a small-cell inner drift chamber were used to record the remaining 
data sample. 

We use a selection procedure similar to that described in Ref. [5j. All charged tracks are required to be consistent 
with originating from the interaction point (IP) ; we impose the requirements dr < 2 cm and | dz | < 4 cm, where dr and dz 
are the impact parameters perpendicular to and along the beam direction with respect to the IP. Particle identification 
requirements are based on CDC, ACC and TOF information [§]. Charged kaon and proton candidates are required 
to be positively identified: the identification efficiencies typically exceed 90%, while misidcntification probabilities are 
less than 10%. No identification requirements are applied for pion candidates, as the pion multiplicity is much higher 
than those of other hadrons. Kg (A ) candidates are reconstructed by combining tt + it~ (p7r~) pairs with an invariant 
mass within 10MeV/c 2 of the nominal Kg (A ) mass. We require the distance between the tracks at the Kg (A ) 
vertex to be less than 1 cm, the transverse flight distance from the IP to be greater than 1 mm and the angle between 
the Kg (A ) momentum direction and its decay path to be smaller than 0.1 rad. Photons are reconstructed in the 
ECL as showers with energies more than 50MeV that are not associated with charged tracks. 

J/ip candidates are reconstructed via the J/ip^£ + £~ (£ = e, /x) decay channel. Two positively identified lepton 
candidates are required to form a common vertex that is less than 1 mm from the IP in the plane perpendicular 
to the beam axis (« 98% efficiency). A partial correction for final state radiation and bremsstrahlung energy loss 
is performed by including the four-momentum of every photon detected within a 50 mrad cone around the electron 
and positron direction in the e + e~ invariant mass calculation. The J/tp signal region is defined by the mass window 
\Mi+£- — mj/^\ <30MeV/c 2 (~2.5c). A mass-constrained fit is then performed for the signal window candidates, 
to improve the center-of-mass (CM) momentum V*ji^ resolution. QED processes are suppressed by requiring the total 

charged multiplicity (A c h) in the event to be greater than 4. In the T(4»S) data J/ip mesons from BB events are 
removed by requiring V*ji^, > 2.0 GeV/c; no requirement on P*j ^ is applied in the off-resonance data sample. 

We also reconstruct charmonia decaying to J/tp- ip 1 candidates are reconstructed via the decay to J/ip> tt + it~ , with 
the ip' signal window defined by \Mj/^ 7! + 7r - — m^i\ < 10 MeV/c 2 (~3tr). Xd and \ci candidates are reconstructed 
using the J/ipl mode; signal windows of ±20MeV/c 2 are chosen around the corresponding nominal masses («2.5c). 
In addition we require cos # 7 < 0, where # 7 is defined as the angle between the photon momentum and the CM system, 
seen from the Xci(2) rest frame. This requirement suppresses the large combinatorial background due to low energy 
photons by more than an order of magnitude, while retaining 50% of the signal, independent of the Xci(2) polarization. 

We use only charged final states for charmed hadron reconstruction to avoid correlated multiple candidates. Can- 
didate D° mesons are reconstructed in the K~ir + , K + K~, Kgir + ir~ and K~ir~ir + 7r + decay modes [10|. We recon- 
struct D + mesons using A _ 7r + 7r + , K~ K + ir + , Kgir + and Kgir + ir + n~ decays; for Df meson reconstruction we use 
the K~K + ir + and Kgir + , and finally A+ baryons are reconstructed via pK~ir + , pKg and A°7r + . A ±15MeV/c 2 
mass window (« 2.5 a) is used throughout, except for the D° — > K~ir~n + iT + and D + — > A'g7r + 7r + 7r~ modes where 
the resolution is better, and the combinatorial background higher: in these cases a ±10MeV/c 2 window is chosen 
(«2.3<r). To study the contribution of combinatorial background under the various charmed hadron peaks, we use 
sidebands selected from a mass window four times as large. 

We generate large Monte Carlo (MC) samples of double charmonium production and of the process e + e~ — > J/ipcc 
with fragmentation to open charm. We also generate a sample of e + e~ — > J/ipqq events for the study of the 
e + e _ — > J/ip X non _ c£ process. In the MC samples the J/ip kinematical characteristics (momentum spectrum and 
angular distributions) are tuned to those measured in the data. As the measured distributions are extracted from 
the data using the MC simulation, the tuning procedure is repeated until the difference between successive iterations 
becomes negligibly small. 

To measure the contribution of cc resonances to the e + e~ — > J/ipcc cross section, we reconstruct all double 
charmonium final states that can result in the presence of a J/ip in the event: J /ip{cc) les , ip'{cc) Tes , and Xci(2)(cc) res , 
where (cc) res is one of the charmonium states below open-charm threshold. If a charmonium state lies above the 



open-charm threshold we assume it will decay predominantly to charmed hadrons; production of a J/ip together 
with charmed hadrons is treated separately below. The process e + e~ — > Y(cc) rcs , where Y is one of the states odd 
under charge conjugation, recently observed in initial state radiation (ISR) studies |12j, can produce J/tp from Y 
decays. However, we are unable to measure this contribution because of the large intrinsic width of the Y states, and 
ignore it. Following the method described in 0, [r| we first reconstruct a (cc) tag = J/^p, tp' , or Xci(2) meson to tag 
the process, and then form the recoil mass 

M rocoil ((cc) tag ) = yj '(Ecu - £ t *ag) 2 " &£, (1) 

where E^ ag and p^ ag are the CM energy and momentum of the reconstructed charmonium, and Ecm is the CM energy. 
The M roco ii((cc)tag) spectra for the data are presented in Fig. [TJ We assume that only charmonium states with a 
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charge conjugation eigenvalue opposite to that of (cc)t a g can appear; two virtual photon annihilation, which can 
produce a pair of charmonium states with the same eigenvalue, was not observed in Ref. [l3| . and is expected to be 
small. 

We fit the four M reco i 1 ((cc) ta g) spectra simultaneously to fix the ipW x c i(2) contributions, which are poorly resolved 
in the M reco u(V>^) spectra. The ratios of the ipv> Xd(2) signal contributions to the M Teco u(ip('*>) and M recoi i(x cl ( 2 )) 
spectra are fixed according to the MC study. The signal line shapes for all the double charmonium final states 
are obtained from MC simulation, with ISR included, and the background is parameterized by a linear function 
(a second order polynomial function in the M leco i\(J / tp) case). Only the region below the open-charm threshold 
(-^recoil < 3.7GeV/c 2 ) is included in the fit. The fitting function for the M recoil ( J/ip) spectrum also includes the 
expected contribution from the ISR process e + e~ — > tp' 7, which is poorly described by the polynomial function; its 
shape and normalization are fixed from the MC simulation. This process was studied for our paper and the 
measured width T ee (ip') was found to be in good agreement with the PDG value 
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FIG. 1: The mass of the system recoiling against the reconstructed a) J/ip, b) tp' , c) Xci and d) Xc2. The curves show the fit 
results described in the text. 



TABLE I: e + e — > (cc)t ag (cc) ICB signal yields (significances) from a simultaneous fit to M reco ii((cc)ta g ) spectra. 



(cc)ros 


J/ip 


(cc)ta g : 

Ip' 




Xcl 




1032 ±62 (19) 


161 ±22 (8.2) 












16 ±5 (3.2) 


9 ±4(2.1) 


XcO 


525 ±54 (9.6) 


75 ±19 (4.3) 






Xcl 


119 ±39 (3.2) 


12 ± 12 






h c 






4±6 


1±5 


Xc2 


99 ±43 (2.1) 


7 ±16 






V'c 


679 ± 63 (10) 


81 ±19 (4.5) 






Ip' 






6±6 


2±5 



The fit results are shown in Fig.[T]by solid curves; the background function and the e + e~ —yip' 7 reflection are shown 
with dashed and dotted curves, respectively. The signal yields and significances for all the studied double charmonium 
processes are listed in Table |U The statistical significance of each process is determined from — 2 \n(£o/£ max ), where 
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TABLE II: e + e — > J/tpH c X signal yields and significances from fits to the M e + e - spectra; the effective reconstruction 
efficiency times branching fraction is also shown. 





D° 


D+ 


Df 


A+ 


N 


1072 ± 108 


715 ± 93 


129 ± 42 


43 ±20 


Significance 


10.1 a 


7.8 a 


3.2 a 


2.2 fj 


Efficiency X B 


0.041 


0.047 


0.022 


0.020 



£ max is the maximum likelihood returned by the fit, and Co is the likelihood with the corresponding contribution 
set to zero. The results for J/ip (cc) ros and ip' (cc) rcs are in good agreement with our previous measurements [Bl. Il3j|. 
There is also evidence for J/ipXci production at the 3.2a level (statistical only). 

Next, we study associated production of a J ftp with charmed hadrons. In the previous paper [j| we determined the 
J/ip cc cross section from measurements of the production rate of a J/ip with associated D° and D* + mesons using 
Lund [l5| model predictions for probabilities of fragmentation cc — ► D°(D* + ). Moreover, to suppress combinatorial 
background from BB events, we applied additional kinematical criteria; the efficiency of these criteria also contributed 
to the model dependence of the result. To eliminate the model dependence in this analysis we use all the ground state 
charmed hadrons: H c = D°, D + , Df and A+, except for Sc ( ' _ ' ) and whose production rates in cc fragmentation 
are expected to be smaller than 1% according to the Lund model. As two charmed hadrons are produced in cc 
fragmentation, the J/ip cc cross section is given by half the sum of the J/tp H c X cross sections. We extract J ftp H C X 
yields in both H c signal and sideband windows, using fits to M e +p- distributions with signal and second order 
polynomial background functions. The J/ip signal shape is obtained from MC simulation, with the small difference in 
the J/tp resolution between the MC and data corrected. The M e + g - spectra are shown for D°, D + , Df and A+ signal 
windows in Figs. [2] a), b), c) and d), respectively; scaled sideband distributions are superimposed. The J/%pH c X 
yields are calculated as the difference between the J/ip yields in the signal window and the (scaled) sidebands. The 
fit results are listed in Table HU We observe a significant excess J/ip signal in the D° and D + signal windows with 
respect to the corresponding sidebands, demonstrating large e + e~ — > J/ip D° (D + ) X cross sections. An excess, with 
low significance, is also seen in e + e~ — > J/ip Df(Af ) X. 

Next, we measure the J/ip momentum spectrum in inclusive production and from the process e + e~ — * J/ip cc. The 
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inclusive J/ip momentum spectrum is obtained by fitting l + i mass distributions in bins of P*j/^ with signal and 
second order polynomial background functions. In the region P*j^ < 2.0 GeV/c only the continuum data is used; the 
J/ip yields are then scaled according to the ratio of luminosities. The ISR processes e + e~ — > tpv) -y contribute to the 
selected sample (with multiplicity greater than four) in the case of fake track reconstruction and/or 7 conversion. 
This contribution is small (~ 2% of the total J / ip rate) , and is subtracted using a MC simulation with ip' and J/ip 
dielectron widths fixed to the PDG values [ljj]. The final yield in each momentum bin, after subtraction of QED 
background, is then corrected for the J/ip reconstruction efficiency and divided by the total luminosity. The result, 
representing the differential cross section, is shown in Fig. [3]a) with open circles. 

We calculate the momentum spectrum of J/ip mesons from all double charmonium processes, including J/ip from 
cascade decays, and note that the final state in e + e~ — > ip> Xci(2) events may contain two J/ip's. We use a MC 
simulation with the contributions of double charmonium processes fixed to the results of the fit to data (Fig. [1] and 
Table [XJ) to obtain this spectrum, shown in Fig. [3]a) with filled circles. The momentum spectrum is peaked near the 
kinematical limit as expected for two-body processes; ISR results in a tail to lower momentum values, and there is an 
additional contribution at £>jy^~3GeV/c due to J/ip's from cascade decays. 

To obtain the J/ip momentum spectrum from the process e + e _ — > J/ip H C X, we measure J/ipH c X yields in bins 
oip*j^. The fits to M^+ e - spectra (Fig. [2] and Table HT]) are repeated in the H c signal and sideband windows for each 
bin, with the J/ip H c X yield defined as the fitted J/ip yield in the H c mass window after subtraction of the scaled 
yield in the H c sidebands. Using the continuum data it is possible to perform such fits below 2 GeV/c, though with 
much larger statistical errors. The yield in each bin is then corrected for the J/ip and H c reconstruction efficiencies, 
using a MC simulation. The sum over all H c weighted by a factor of 0.5 is plotted in Fig[3]a) with filled squares and 
represents the J/ip momentum spectrum from the process e + e~ — > J/ipcc, where the cc pair fragments into charmed 
hadrons. The sum of this distribution and that from double charmonium production represents the J/tp momentum 
spectrum from the process e + e~ — > J/ipcc; it is shown in Fig.[3]b) by the open squares. The difference between this 
and the inclusive J/ip spectrum is thus the spectrum from e + e _ — > J/ipX non _ cB events, where the system recoiling 
against the J/ip is not produced via a cc pair (shown by the filled triangles in Fig. [3]b), to which the color-singlet 
e + e _ — > J/ipgg and color-octet e + e~ — > J/ipg processes contribute. 




FIG. 3: J/ip momentum spectra: a) inclusive (open circles), from e + e~ — > J/ipH c X (filled squares) and from double char- 
monium production (filled circles); b) the sum of all e + e~ — » J/ipcc processes (open squares), from the e + e~ — > J/ip X non . C c 
processes (filled triangles). The results of fits to the Peterson function are shown in a) for the inclusive spectrum (solid curve); 
in b) for the processes e + e~ — > J/ipcc (solid curve), and e + e~ — > J/ip X non . C c (dashed curve). 
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TABLE III: Cross sections for the processes e + e — ► J/ip X, J/ipcc and J /if) X non . C c ([pb]), and characteristics of the J/ip 
spectra (ep e ter, Qhci and a pr od); 'X / n dof values for the corresponding fits are listed in parentheses. 





J/ipX 


J/ipcc 


J J Ip -Xnon-ee 


a 


1.17 ±0.02 


0.74 ± 0.08 


0.43 ± 0.09 


O" Peter 


1.19 ± 0.01 


0.73 ±0.05 


0.48 ±0.07 


f Peter 


0.16 ±0.01 (8.9) 


0.10 ±0.02 (0.6) 


0.32±g;il (1.6) 


Qhel 
Q^prod 


0.03±0.03 (0.6) 
0.69±0.05 (3.3) 


-0.19±lj;i (1.0) 
-0.26t° ; 2 22 (0.5) 


0.4ll u ;S (1.2) 
5.2±l;i (0.3) 



The J/ip efficiency corrected momentum spectra, shown in Fig. [3] for the processes e + e~ — > J/tpX, J/ipcc and 
J /if) X non - C c, are then used to calculate the respective cross sections, after performing a sum over all momentum 
bins. The results are presented in Table IIIIl The statistical errors are dominated by the momentum interval 
Pj/ip ^ 2.0GeV/c, where only the small continuum data sample is used. To characterize the hardness of the mo- 
mentum spectrum, we perform fits using the Peterson function [16J; the parameters ep e ter for the e + e~ — > J ftp cc and 
J j 'ip X non _ cS processes are listed in Table [TTT1 For completeness, the resulting cross sections <7p e ter are also shown: they 
are consistent with the directly calculated values, with statistical errors reduced by a factor of 1.5, as the fit effectively 
extrapolates the high- momentum results into the low- momentum region. Such results are model-dependent, and we 
rely instead on the directly calculated values a for the cross section. 



We note that unlike our first paper 17j no correction for the iV c h requirement is applied for any of the process 



studied. For e + e~ — * J/tp X non _ cS such corrections are only possible by relying on a model. However, for the process 
e + e~ — * J/ipcc, the efficiency of the N c h > 4 requirement is more than 99% if the cc pair fragments into charmed 
hadrons, as their decays lead to a large multiplicity in the event. For double charmonium production the efficiency is 
70% according to the model used in the MC generator, and varies by ±20% with different charmonium decay models. 
As double charmonium represents only ~ 10% of the total e + e~ — > J/ipcc cross section, the resulting correction is 
small, and included in the systematic error. 

We also perform an angular analysis for the e + e~ — * J/ipcc and e + e~ — * J/ipX non . c5 processes. This provides 
important information on the production mechanisms, and allows the efficiency calculation to be improved: the J/ip 
reconstruction efficiency depends on both the production angle (# pr od, the angle between the J/ip momentum and the 
beam axis in the CM frame) and the helicity angle (#hci, the angle between the £ + from J/ip decay and the CM, seen 
from the J/ip rest frame). The MC simulation is adjusted to match the measured distributions. 

Angular distributions are obtained from fitted yields in bins of |cos(9 pro d| and |cos(9h c i|, with an appropriate effi- 
ciency correction performed bin- by-bin, for inclusive J/ip, J /ip from double charmonium production, and J/ip from 
e + e~ — > J/ipH c X. The results are shown in Fig.[U The inclusive J/ip distributions (open circles) are obtained from 
J/ip yields. Those for double charmonium production are obtained from fits to the four M reco ii((cc)t a g) distribu- 
tions, as for Fig. Q] above. Distributions for e + e~ — * J/tpH c X are obtained from fitted J/ip yields in appropriate H c 
mass windows, after subtraction of yields in the H c sidebands. The distributions for e + e~ — * J/ipcc (open squares) 
are calculated as the sum of the corresponding distribution for double charmonium production (with weight 1.0) 
and e + e _ — > J/ipH c X (with weight 0.5). Distributions for the e + e~ — > J/ip X non . c5 process (filled triangles) are 
determined from the difference between e + e~ — > J/ipX inclusive and J/ipcc distributions in each bin. 

We fit the helicity angle distribution with a function ~ (l + c%el cos 2 (#i lc i)). While the production angle distributions 
are also fitted with a function ~ (1 ± a pro d cos 2 (# pro d)), we note that these distributions can differ from 1 + a cos 2 9 
due to ISR or the contribution of the e + e~ — ► 7*7* — > J/ipX process to the J/ip X non _ c5 final state. The fits yield 
the parameters cehei and a pro d listed in Table IIIIl the fit results are shown in Fig. 2] 

The systematic errors on the production cross sections for both e + e~ — > J/ipcc and J/ip X non _ ce processes are 
summarized in Table IIVI In the double charmonium production study, systematic errors due to J/ip yield fitting are 
determined as in our previous papers 0, [T3] ] ; we also perform variant fits including final states with two charmonia 
with the same charge conjugation eigenvalue. In the study of associated production, we consider changes in J/ip H c X 
yields under variation of the fitting procedure (a 2-dimensional fit to (M(J/ip), M(H C )), a fit to the M(J/ip) distri- 
bution in bins of M(H C ), and to the M{H C ) in bins of M(J/ip)), as well as variation of the signal and background 
parameterizations, the fit ranges, and the binning. The uncertainty in H c reconstruction efficiencies due to the un- 
known kinematics of cc fragmentation into charmed hadrons is small, due to the weak dependence of reconstruction 
efficiency on H c momentum, and is included in the total systematic error. 
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FIG. 4: Angular distributions (jcos#hei| in a), |cos6lp ro d| in b)) for inclusive e + e~ — * J/ipX (open circles), e + e~ — > J/ipcc 
(open squares), and e + e~ — > J/-ip X non . C c processes (filled triangles). The results of the fits described in the text are shown 
with the dash-dotted, solid, and dashed curves respectively. 



TABLE IV: Summary of the systematic errors on the cross sections shown, in percent. 



Source 


J/ipX 


J/4>cc 


J j ijj X n on-cc 


Fitting procedure 


±3 


±5 


±9 


Selection 




±3 


±5 


Angular distributions 


±4 


±6 


±10 


iVch requirement 




+5 
-0 




ISR 




+4 
-2 


+4 
-7 


Track reconstruction 


±2 


±5 


±8 


Identification 


±2 


±4 


±7 


B(J/i>), B(H C ) 


±1 


±3 


±3 


Total 


±6 


+ 12 
-11 


±20 



When the integral J/ip production and helicity angle distributions in the MC simulation are tuned to those in the 
data, their correlations are not taken into account. We assume the most conservative correlations, resulting in the 
largest deviation of the J/ip reconstruction efficiencies that reproduce the integral distributions. The resulting differ- 
ence in efficiency is the largest contribution to the systematic error. Other contributions come from the uncertainty in 
the track and Kg (A ) reconstruction efficiencies; from lepton, kaon and proton identification; and from uncertainties 
in absolute H c branching fractions. 

In summary, we have measured the cross sections for the processes e + e _ — > J/ipX, J/ipcc and J/ipX non _ cS to 
be (1.17 ± 0.02 ± 0.07) pb, (0.74 ± 0.08 torn) P b and (°- 43 ± °- 09 ± °- 09 ) P b ' respectively. We therefore conclude that 
e + e~ — > J/ipcc is the dominant mechanism for J/ip production in e + e" annihilation, contrary to the expectation from 
NRQCD. Moreover, this cross section exceeds the perturbative QCD prediction er(e + e~ — > cccc) «0.3pb [HI], which 
includes the case of fragmentation into four charmed hadrons, rather than J/tpcc. The cross section for J/ip X non _ c£ , 
which can proceed via e + e~ — * J/ip gg or J/ipg, as well as e + e~ — > J/ipj* diagrams, is of the same order as that for 
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J/ip cc. We have measured the J/tjj momentum spectrum and the production and helicity angle distributions from 
all three processes. For the e + e~ — > J/tp X noa . C c process, the J/tp momentum spectrum is significantly softer than 
that for e + e _ — * J/ipcc, and the production angle distribution peaks along the beam axis. We note that all the 
measured cross sections are full (rather than Born) cross sections and include contributions from cascade J/V', and 
that model-dependent corrections for the charged track multiplicity requirement have not been performed. 
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